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VENTILATOR PATIENT SYNCHRONIZATION 

This application claims the priority filing date of U.S. provisional patent 
application serial no. 60/292,983 filed on May 23. 2001 . 

5 FIELD OF THE INVENTION 

This invention relates to methods for synchronizing ventilators with patient 
respiration. More specifically it relates to methods for estimating the phase in patient 
respiration using patient flow measurements so that the respiratory support of a 
ventilator can match the phase of the respiratory cycle of the patient. 

10 BACKGROUND OF THE INVENTION 

Mechanical ventilators assisting spontaneously breathing patients strive to 
synchronize their perfomiance with the patient's efforts. To do this, ventilators 
typically measure one or more of pressure, volume, flow and time and compare that 
measure with predetermined thresholds. Some ventilators use respiratory bands 

15 around the chest and abdomen of the patient to detect patient effort:. The ventilator 
may then adjust the pressure, volume or flow of air being delivered to the patient in 
accordance with measure of the patient's efforts. For example, a flow-triggered 
pressure controlled device may deliver air at one fixed pressure to a patient until the 
flow crosses a threshold level, whereupon the pressure is changed to another fixed 

20 pressure. Depending on their conditions, different patients may experience different 
levels of discomfort depending upon how quickly and accurately the ventilator tracks 
the patients' efforts. Simple threshold tests may fall when breaths are irregular, for 
example, during the presence of coughs, sighs and snores. An improved method 
and apparatus for ventilator patient synchronization is described In Patent 

25 Cooperation Treaty Applfcation PCT/AU97/00631 with publicatfon number WO 
98/12965 (Berlhon-Jones) where the phase of the patient's respiratory cycle is 
determined from flow data using fuzzy logic. The specification is hereby included by 
cross-reference. 

Typical apparatus includes a sen^o 3 controlled blower, comprised of a motor 
30 2 and an impeller 1 connected to a patient interface 5 via an air delivery conduit 6, as 
shown in Rg. 1. The controller 4 is typically a computer, a processor including 
memory, or a programmable circuit. One example of patient interface is a nasal 
mask, others include nose and mouth masks, full face masks and nasal pillows. The 
pressure in the mask may be measured by a transducer 1 1 having direct contact with 
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the mask, or alternatively, the transducer may be physically situated in the blower 
main housing and may estimate the mask pressure using correlations. Flow 
transducers 1 0 or other means for measuring flow may also be situated in the mask 
or in the blower main housing. There are various displays 8 and switches 7 on the 
5 blower housing. There is an interface 1 5 to enable the apparatus to communicate 
with other devices. Some apparatus include a fixed speed blower whose output is 
controllably variably vented to atmosphere providing a controlled variable pressure to 
the patient. 

"Expert" systems are known to be used for assisting with medical diagnosis. 

10 Such expert systems are lypkally said to comprise two parts, a "knowledge base" 
and an "infterencing engine." The knowledge base comprises the set of "expert" 
information about the system which is used to guide interpretation of the data which 
has been obsen/ed. Sophisticated expert systems may include hundreds, or 
thousands of pieces of information in the knowledge base. The fuzzy membership 

15 njles and weights of Berlhon-Jones may be interpreted as the knowledge base. The 
inferencing engine is the mechanism which combines the knowledge base with the 
experimental evidence to reach the conclusion. Several different inferencing engines 
are known, such as those based on fuzzy logic, rule based reasoning and Bayesian 
likelihoods. 

20 Bayes' theorem^ quantifies the intuitively appealing proposition that prior 

knowledge should influence interpretation of experimental observations. One form of 
Bayes' theorem is: 

25 where L Is a likelihood or probability function. Thus, UH|F) Is the likelihood of an 
hypothesis being true, given observatton F. L(H) is the likelihood of the hypothesis 
being true, and UF|H) Is the likeBhood of the observation given the hypothesis being 
tme. 

For example, if a physician has observed a particular symptom in a patient, in 
30 deckling whether the patient has a particular disease, the physician draws upon the 
prior evidence of the likelihood that the patient has the particular disease. Several 
independent observattons may be used In conjunction with prior likelihoods to 



wo 02/094358 



PCT/AU02/00639 



3 



deleimine the likelihood that an hypothesis is true. The decision may be taken to be 

the most likely hypothesis. 

BRIEF DESCRIPTION OF THE INVENTION 

The invention is a method and apparatus for detemr^lning phase with 

5 probability functions. The method Involves partitioning the respiratory cycle into 
discrete phase states. The states will include inspiration and expiratfon and 
preferably include a number of additional states within inspiiaton and expiration. The 
probabilities of each phase slate are then calculated using probability functions. In 
the preferred embodiment, the cateulaBon Is a function of an "observed prDbability" 

10 determinatfon H^Hd and a Drior probability" detemiination L(Hi). The calculated 
probabilities L(H,|Fi) are then compared to detemtiine the actual phase in the patienf s 
respiratory cycle that the patient is experiencing. 

Prior probabilities of each particular state may be detennined by a function of 
relative duration, that is. the ratio of each state's duration to the duration of an entire 

15 respiratory cycle. The invention contemplates discrete phase states that may be 
equivalent in duration or of unequal durations. Such states of different duration may 
be derived by partitioning the respiratory cycle into sections in relation to a peak ftow 
or differing rates of change of flow. Moreover, cateulated probabilities my be adjusted 
based upon state successton. e.g.. increasing the likelihood of the slate that follows 

20 consecutively from the cunently determined state. 

The obsen/ed probability determination evaluates the likelihood of each state 
given an obsen^ed respiratory characteristfc. In the preferred embodiment of the 
invention, respiratory flow Is observed or measured from the patient In the 
evaluatton, the measure Is applied to a probaWllty function that Is based upon the 

25 predetermined likelihoods between that measure and each potential phase state, for 
example, the likelihood that a particular range of flow would be associated with a 
particular state. 

•nie method may also be used to dynamically improve reliability of phase 
detection by learning the respiratory patterns or characteristics of the patient. By 
30 utilizing previously calculated state probabilities to adjust prior probabilities. 

subsequently cakjulated state probabilities provide greater patient specific accuracy 
in phase detemiination. Moreover, previously recorded breath data can be used to 



' Armitaae & Berry (1994) Statistical Methods in Medical Researdi. 3rd Edition. p72. 
bSSs^ Ud. Oifo«l. united Kingdom ISBN 0^-03695-8 
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modify the partitioning of phase states or revise the relationship between observed 
respiratory characteristics and the likelihoods for each state given the observation. 

BRIEF DESCRIPTION OF DRAWINGS 

Rg. 1 depicts a prior art ventilator apparatus suitable for use with the invention; 
5 Fig. 2 is a flow chart of steps in a method in accordance with the invention; 

Rg. 3 illustrates phase classification based on equivalent stale duration; 

Rg. 4 Illustrates phase classification based upon peak flow; 

Rg. 5 shows phase classification with a set of obsen/ations useful for determining 

probabilities of the identified states; 
10 Rg. 6 shows flow versus time data for a person over a period of 1 minute; 

Rg. 7 shows flow versus time data for a person over a period of 10 seconds. 

DETAILED DESCRIPTION OF THE INVENTION 

The invention provides a mettiod and apparatus for improved ventilator 
patient synchronization In which the most likely phase is detemiined on the basis of 

15 flow and rate of change of flow infomiaflon using Bayes' tiieorem or some similar 
formula for assessing conditional probability. In an identifying step 20 as shown in 
Rg. 2. according to the invention, the respiratory cycle Is divided or partitioned into a 
number of distinct states, for example, early inspiration, mid-inspiration, late 
inspiration, earty expiration, mid-expiration, late expiration, pause, sigh and cough. 

20 When the metiiod and apparatus according to the invention are being used, flow and 
pressure obsewations are made continuously and at each time point in a measuring 
step 22. In a calculating step 24, probabilities or likelihoods are detemiined. The 
likelihood of tiie phase corresponding to each of tiie states is detemnlned preferably 
using Bayes' theorem as shown in equation 1 , supra. For example, in a system that 

25 distinguishes between six states, six calculations for each distinct state using 

equation 1 would be made. This evaluation may include an assessment of tine prior 
likelihood of each of the identified states L(H,). However, these may be detennined 
before use and/or adjusted during use. Ukelihoods of tiie distinct respiratory states 
occun-ing given flow and rate of change of flow observations L(Hi|Fj) are also 

30 detemiined. After tiiese evaluations, in a detemnining step 26. the patient's phase is 
detemiined to be the phase which con-esponds to the most likely state. e.g. the phase 
' cafculation with tiie highest probability or likelihood. This Infomiation may then be 
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used by the apparatus to adjust the flow, volume or pressure of air being delivered to 
the patient in order to improve ventilator synchronization. 

The prior likelihood or prior probability of each state L(H,) depends upon the 
definition of the state and its relative duration. Such a determination may be 
5 accomplished by a ratio of the duration of the state to the duration of a respiratory 
cycle. Each state may not necessarily have the same duration. A simple system 
might divide the respiratory cycle Into two slates: Inspiration and Bcpiralion. 
Alternatively, a greater number of states within inspiration or expiration may be 
classified by the system. If the respiratory cycle ted to sinusoidal flow vs. time 
10 curves, then it may be appropriate to divide the cycle into equal temporal subsections 
and assign each subsection as a state in which case, each slate may be equally 
likely a priori. For example, if the duration of a hypothetical respiratoiy cycle was 6 
seconds and the system distinguishes six slates, then each state may be one second 
in duration, the first second being temned "eariy inspiration", tiie next second "mid- 
15 inspiration" and so on. In this model, ttie different states may have the same prior 
likelihoods. TTius, as illustrated in Fig. 3 the likelihood of each state L(H^ is 1/6 as 
follows: 

L(HeailyJnspiralion) = 1/6 
20 L(HmkJJnspi™ilon) = 1/6 

UHiateJnspIratton) = 1/6 
L(Heaily_explialion) = 1/6 
L(Hirt(Lo)(plratlon) =1/6 
L(HiaiB_ex|)lmaan) = 1/^ 

Increasing ttie number of states increases the potential accuracy of a system 
to track changes in respiratory state, however, this comes at an increasing 
computational cost. It may be computationally more efficient to have states of 
different duration. In one model of respiration, ttie states may be chosen on ttie basis 

30 of magnihide of flows, for example earty flows from 0 up to 1/eth of tine peak flow may 
be tenned "eariy inspiration", in such a system, for example, as shown in Fig. 4, 
again based upon a respiration cycle for a typical patient flow L(HeariyjnspiraHon) might 
be 4/48 where 4 approximates tiie duration of ttie state and 48 approximates the 
duration of tiie entire cycle, viewed as ttie intercept of ttie peak flow range and ttie 

35 flow cycle. Similarly, a mid inspirational state detemiined from 5/6 to 6/6 of tt.e peak 
flow might have a likelihood L(H^j„s,,aton) equal to 1 1/48. Fig. 4 also shows a mid 
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expiration state with L{H^.^ equal to 9/48. The remaining states may be 
determined in similar fashion. 

In another model of respiration, the rate of change of flow may be chosen as 
the basis for dividing the respiratory cycle. Thus, an early inspiration state may be 
5 characterized by flow where the rate of change Is in a range appropriate to classify 
the rate of change of respiration in its eariy stages. Using for example a graph, the 
duration of a discrete portion of the respiratory cycle in that particular rate of change 
range divided by the duration of the entire cycle can then yield the prior probability of 
the particular state L(Hi). An advantage of dividing the respiratory cyde on the basis 
1 0 of rate of change of flow is that those parts of the respiratory cycle where flow 

changesmorerapidiymayhavemorerespiratofyslatesasslgnedtothem. Inthese 
latter two models, the different slates may have different priorlikellhoods. 

In one torn of the invention, the number of slates and their definitions may be 
adaptively refined depending upon the breath data. For example, where the flow rate 
15 is changing rapidly, the apparatus according to the invention can detennine that it is 
appropriate to have more states to define that region. 

Obsen«»tlon of breath data of patients over several nights leads to the 
conclusion thatacough is aprfon less likely to occur than a breath, since patients 
breath more often than they cough. It may also be observed that a deep breath is 
20 more likely following a number of shallow breaths. Patients are known to pass 

through a number of different stages during sleep, such as REM sleep. The nature of 
their breathing can change depending on the different stage of sleep and hence 
different types of breath are a priori more likely depending on the stage of sleep. To 
the extent to which a particular patient exhibits characteristic breaths, the system has 
25 the capability to learn these characteristics and adjust the prior likelihoods of different 
states to suit the patient as herein described. 

A discrete set of /observations for purposes of determining an obsenred 
probability function UFjH.) may be defined for a flow signal. SimilarV. a discrete set 
of observations of the derivative of flow may be defined. The function may be 
30 predefinedorchangedasaresuKofchanglngpatientraspiratorycharacteristics. The 

function may be defined by dividing the range of values for flow or derivative of flow 
into discrete sub-ranges and then comparing those sub-ranges to each state and 
assessing probability based upon a typical or pieviously obsen/ed respiratory flow 
ounce's relationship with that sub-range. For example, the flow signal may range 
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from -200 l/min to +200 l/min. This may be divided into the following sub-ranges: 
200 to -50, -50 to 0. 0 to 50. 50 to 200 l/mln. Noise in the flow signal defines the 
smallest division into which flow signals may be divided (i.e., there is no point in 
dividing the total range of flow signals into 1 l/min divisions if the noise level is 2 
l/min). The following table, labeled Observation Table 1 , gives an example of the 
likelihoods of the obsen/ations being made, given the various hypotheses being true 
for a system that distinguishes six states. In this example, there is a low likelihood 
that a large positive flow signal will be observed during the mid e)9)iration phase. 
OBSERVATION TABLE 1 



L(F|H) 


Early 

Inspiration 


Mid 

Inspiration 


Ute 

Inspiration 


Early 
Expiration 


Mid 

Expiration 


Late 

Expiration 


200 to 100 

(Large 

positive) 


low 


high 


low 


low 


low 


low 


50 to 100 
(Medium 
positive) 


medium 


medium 


medium 


low 


low 


low 


0to50 
(Small 
positive) 


high 


low 


high 


medium 


low 


medium 


-50 too 

(Small 

negative) 


medium 


low 


medium 


high 


low 


high 


-100 to -50 

(Medium 

negative 


low 


low 


low 


medium 


medium 


medium 


-200 to - 
100 (Large 
negative) 


low 


low 


low 


low 


high 


low 



10 



While Obsen/ation Table 1 includes the terms "low", "medium" and "high" to 
describe the likelihoods or probabilities, in actual calculations these temfis vwll be 
quantified in the ranges of 0-33 percent for low, 33-66 percent for medium and 66 - 
100 percent for high. Moreover, these probabilities need not be limited to three 
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levels. Generally, a change in the boundaries of the states leads to a change in the 
likelihood of each obsen/ation given the hypothesis that the patient is in a given 
respiratory state is tme. Alternatively, the probabilities may be determined from 
defined continuous likelihood density functions. 

5 Rg. 5 shows sample breath data that is subdivided into 9 states (sho\wn with 

Roman numerals). Note that states "i", °ii" and V are relatively brief compared to 
states "iii" and "vi". Hence, everything else being equal, states "lii" and "vl" have 
relatively greater prior likelihoods than states "i". "B" and V. The figure also depicts a 
set of flow observations (large positive to large negative) as the obsen^attons relate to 

10 the classiflcaHon of phase states. The observed probabilities UFi\Ht) may be 

determined from this graph. For example, the Intercept or observation of a flow in the 
range of 80-100 (large positive) would suggest a high probability for state ill but tow 
for states i. v, vi. vil, viii, ix. States ii and iv may have low to medium likeNhoods as 
these cases are at the boundary of the large positive observation. 

15 In the following hypothetical example, a calculation of the phase probability or 

likelihood is demonstrated In accordance with the invention utilizing the simplified six 
dasslfteation system of Fig, 3 in detemiining L(Hi) and Observation Table 1 for 
determining L(Fi|HO. The results are as follows: 



20 Observed fact F: 

Measured flow from ftow signal is 150 

So Fj = Fiaigej)osMve 

25 Ukelihood of hypotheses UHiY . 

UH«.,u™pi«fion) = 1/6 = 0.167 
L(H,rt<UM«iio,>) = '>/6 = 0.167 

UHiateJnspImtkx.) = 1/6 = 0.167 
30 L(Heaily_explratIon) = 1/6 = 0.1 67 

L(H,rf<Loxpira««.) =1/6 = 0.167 
L(Hbte_«.p.™«o„) =1/6 = 0.167 

Likelihood of observation given hypothesis L(Hi|Fj) 

KJFm^j^e |Hearty_tasplra«on) = 0.25 (aSSUfTling 0 25 = lOW) 

L(F,a^^. IH^^nspimHoo) = 0-75 (assuming 0.75 = high) 

L(Flarge_positlve [Hiatejnspi ration) = 0-25 
L(F|aigB_posltive |Hearty_e;(piratk)n) = 0-25 
40 L(Fbigej)Osaive |HmkLe)(plfBBon) = 0-25 

L(Fiaig8_posiav© lH|atB„«xplraHon) = 0-25 
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Bayes' Theorem Calculations: 

0.25x0.167 

L(HeartyJnsr^Ban|Flarga_pc>slllva) 0,25 x 0.167 + 5 X (0.1 67x0.25) 



0.125 



0.375 



0.75x0.167 

MH„wjAspj«itonlFte,ge.posiBi.o) =0,25 X 0- 1 67+ 5 X (0.1 67x0.25) 



10 ' MHteteJnspteUFterge^Jfiv^) = 0.25x0.167 = 

0-125 0.25x0.167+5x(0.167x0.25) 

0.25x0.167 ^ 

0 125 MHearty_expf«tlan|F=largejx«^ = 0.25 XO. 1 67 + 5 X(0. 167x0.25) 

15 0.25x0.167 

0 125 '■^"'""-^'^^'^'^■^^ "0.25x0.167+5x( 

0.25 x0.167 

20 0.125 '''''-^'*"^'^0.25x0.167 

By comparing the calculated probabilities L(H,|Fj) to find the most probable phase, it 
Is apparent that in this simplified example the phase is early inspiration. With this 
determination, a system providing ventilatory support may automatically adjust the 

25 support to accommodate the patienf s phase, for example, by adjusting a pressure 
gain associated with the particular state of the detennined phase. Moreover, the 
probability for the phase L(H,) may be modified after calculations to improve 
subsequent determinations by adjusting L(H|) to equal or appn^ach the detennined 
phase likelihoods L(Hi|Fj) by a function of L(Hi|Fj). Moreover, the likelihoods of each 

30 of the obsen/ations L(Fj|Hi) may be updated based upon recorded breath data 

representing a prior respiratory cycle or some average of the prior respiratory cycles. 

In accordance with the invention, a processor of a ventilator apparatus like 
that Illustrated in Fig. 1 is programmed with an algorithm that accomplishes the 
following prefened method:: 

35 I Repeat the following with a first time constant 

(a) Update the prior likelihood for each state L(H|) 

(b) Update the likelihood of each obsen/ation, given each state L(F,|H,) 
II. Repeat the following wrth a second time constant 

(a) Measure flow, detennine rate of change of flow 
40 (b) Calculate ttie likelihood L(H,|F^ for each state using Bayes' ttieorem. 

(c) Phase is determined to be the most likely state. 
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The first time constant may be in the order of several typical breaths, for 
example 10 to 15s. The second time constant is shorter, for example, 1/20th of the 
duration of a typical breath, or several hundred milliseconds. 

In one fonri, the apparatus according to the invention provides a training 

5 period when the patient is allowed to wear the mask in a relaxed environment The 
flow and pressure in the mask are monitored and the device has the opportunity to 
learn some of the characterislics of the patient's respiratory patterns. During this 
period adjustments to the functions for likelihoods P(H,) and P(Fj|H,) for each state of 
the system can be made using L(H,|F,) and obsen/ed breath date and/or the 

10 partitioning of phase states may change. The respiratory characteristics of the 

patient will change when the patient falls asleep, however, the apparatus can monitor 
the transition to the sleep stete and estimate probable parameters for the sleep stete. 

While patient respiratory flow can be modeled as being generally sinusoklal, 
in rBality, the data are more complicated. For example, figures 6 and 7 show flow 

15 date from a person over periods of 1 minute and 10 seconds respectively. The date 
can be additionally complicated with the presence of noise, artefact, drift, sighs, 
coughs and other events. During the presence of flows of small absolute value, the 
effect of heart beat (the so-called "cardiogenic effecf ) on flow data can be identified. 
These complications, which may be termed "non-Idealities", mean fliat a simple 

20 automatic system may incon-ectly decide phase. 

For example, with reference to Fig. 7, a system based purely on tiireshold 
values might conclude that the patient had moved to the inspiratory phase at 
approximately 2 seconds, and then a short while later, for example a hundred 
milliseconds, gone to the expiratory phase, and then switched back to the inspiratory 

25 phase once again. In a simplrfled automatic system, such a breath may complicate 
cataulations of the length of a respiratory cyde and the volume of air inhaled or 
exhaled. Hence there is a need for a method and apparatus which can accurately 
determine the phase of a respiratory cycle of a patient from complex flow date. There 
is a need for a method and apparatus which is tess likely to be confused by non-ideal 

30 flow when detemiining the phase of the patienL 

One approach to solving tiie problem of dealing with non-idealifies in flow date 
is to smooth the date, lor example by tow-pass filtering the date in the measuring step 
22. A low pass filter may be characterized by the time constent of the fitter. A small 
time constent may mean that the amount of smoothing is small, but the system can 
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respond more quickly. A larger time constant may mean that the data is smoother, 
but that the system does not respond as quickly. A very large time constant - 
essentially a long term average - may remove all phase information. Hence there is a 
need to choose an appropriate time constant if the data are to be smoothed. 
5 Furthemiore, there may be an advantage in a system that learns from the flow data of 
the patient what an appropriate time constant may be. Such a system can 
dynamically update the time constant if necessary. 

Detennining the derivative of noisy data can be parficularfy problematic. An 
algorithm for calculating a derivative may give very different results depending on tiie 

10 period over which the data is analyzed. For example witti reference to Fig. 7, an 
algorithm which calculated the derivative of data over ttie region from 2 seconds to 
2.1 seconds may give a negative value, while an algorithm which calculated Vne 
derivative over tiie region from 2 seconds to 2.5 seconds may give a near zero value. 
A turning point is a part of a curve where the slope of ttie cun/e changes between 

15 positive and negative, ttiat is, ttie second derivative is zero. Similarly, the data at ttie 
end of a respiratory cycle (that is, after exhalation and before inhalation) can be 
particularly non-ideal because of the presence of cardiogenic flow. As a result, the 
cateulations for detennining Vne derivative of flow at ttie end of a respiratory cycle can 
be problematic and error prone. 

20 To overcome this problem, an algorithm for detennining phase in accordance 

with the invention may first detennine whetiier ttie patient is likely to be near ttie end 
of a respiratory cycle (for example, if ttie absolute value of a flow datum is small, 
hence distinguishing from peak inspiration and peak exhalation). If ttie patient is 
likely to be near ttie end of a respiratory cycle, ttien a different time constant may be 

25 applied to ttie smoothing of data for tiie purposes of calculating a derivative. In a 
preferred embodiment, ttie time constant will be longer In ttie region near ttie end of 
ttie respiratory flow curve ttian at ottier regions, such as eariy to late inspiration and 
eariy to late expiration. Furthermore, in an algortUim In accordance witti the 
invention, ttie effect of ttie derivative in detennining ttie phase of tiie patients 

30 respiratory cyde is given less weight in ttie region where derivative calculations are 
' likely to be error prone, such as one or more of ttie following: ttie end of ttie 
respiratory cycle, peak inspiration, and peak expiration. Furttiermore. in an algorithm 
in accordance wtth an embodiment of the invention, an a priori likelihood function for 
the derivative of flow in ttie end of ttie respiratory cycle has a maximum at a low or 
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zero value of flow and a minimum at large values of the derivative. In other words, in 
an algorithm in accordance with an embodiment of the invention, in this region, a low 
value of the derivative has a higher likelihood than a high value of the derivative. 

In a preferred embodiment of the invention, 12 bit flow data are collected at 50 

5 Hz. The data are low pass filtered with a time constant of 160ms. 

As previously described, in a system in which phase Is divided Into several 
states, an algorithm in accordance with an embodiment of the invention may 
detemiine that the patient is in one of those states, for example the state with the 
largest likelihood or probability. Alternatively the algorithm may detemnine that phase 

10 is part way between two adjacent states with high likelihoods. The algorithm may 
determine that the phase is closer to the one of the two states with the highest 
likelihood, for example, by weighting in accordance with the relative size of the 
likelihoods. Hence if one slate is 95% likely and another Is 70% likely, the phase will 
be determined to be closer to the state which Is 95% likely. In this regard, phase may 

15 be detenmined as a continuous variable. For example, phase may vary from 0 to 1 
depending on the state determined or depending on the proximity between two most 
likely states. Alternatively or additionally, phase may be represented using polar 
coordinates and vary from 0 to 360 degrees, or 0 to 2 pi radians. 

To illustrate such a detemfiination, phase may be calculated as a continuous 

20 variable ^ . One way to estimate this from a series of posterior probabilities i} of N 
states S, is to associate each state with a standard phase and weighting factor 
JVj then calculate <^ as 

Another approach is to use the formulation of Bayes theorem in temis of contnuous 
25 random variat)ies, which yields the following posterior density function (assuming 
here that the Information available alxjut the patient is the instantaneous respiratory 
flow q), 

and the estimate of ^ is that value which maximizes the posterior density given the 
30 obsen^ed value q of patient flow (this may be a locaJ or global maximum). In 

practice, the prior and conditional densities which are used to calculate the posterior 
density may be estimated from discrete data, possibly by using continuous 
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interpolating functions sucli as piecewise linear or cubic spline functions. The 
posterior density function may then be calculated (possibly piecewise) analytically or 
via numerical methods, yielding a posterior density function which may thus be a 
continuous function of ^ and q. If the posterior density is a continuous function, it 
5 follows that since flow ? is a continuous function of time tiien tiie estimated phase <|) 

is a continuous function of time. 

Such a continuous phase variable may tiien be used to adjust ventilator 
support pressure. For example, if the pressure delivered to tiie patient Is based on a 
continuous function n of phase, 

10 p=p,.«.+aH^) 

where and A are constants or slowly varying quantities, then an important 

consequence of the phase being a continuous function of time is tiiat tiie pressure is 
a continuous function of time, which may be desirable for example in order to avoid 
tiie discomfort which may be associated with sudden changes in pressure delivered 

15 to the ainway. 

Of course, as the number of discrete states increases, the representation of 
phase may be regarded as a continuous, ratiier ttian discrete function. This may 
alleviate a need for the additional calculation of a continuous phase variable from the 
detemiined phase states. In this regard, although six and nine states have been 

20 identified in prior embodiments, by way of example, phase may be divided into as 
many as 128 or 256 states. 

In an algorithm in accordance witti an embodiment of the invention, phase 
determinations may be enhanced by other phase related likelihoods. For example, 
consecutive phase slates are more likely than non-consecutive phases. Thus, a 

25 phase detennination may be based in part upon tiie likelihood tiiat ttie current state is 
the next consecutive state from a previously detemiined state. To illustrate such a 
process, suppose that it has been determined that ttie current phase state is early 
inspiration. If the next phase state is mid-inspiration, tiien mid-inspiration is a priori 
mor© likely than a non-oonsecutive state, such as late inspiration or an expiratory 

30 state. Hence as phase progresses tiirough tiie cycle, a priori likelihoods are adjusted 
in accordance witti their consecutive status, for example by increasing the calculated 
probability of the next consecutive, state. This is not to say that an algorithm in 
accordance with an embodiment of ttie invention cannot make a determination tiiat a 
non-consecutive state will follow tiie current state, but simply that it is less likely to 
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make such a determination. Hence In an algorithm in accordance with an 
embodiment of the invention, phase "cycling" - where the algorithm flips back and 
forward between two states - is less likely. 

Although the invention has been described with reference to partfcular 
embodiments, it is to be understood that these embodiments are merely illustrative of 
the application of the principles of the Invention. Thus. It Is to be understood that 
numerous modifications may be made in the illustrative embodiments of the Invention 
and other arrangements may be devised without departing from the spirit and scope 
of the invention. 
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Claims 

1 . A method for determining a subjecTs respiratory phase during ventilatory 
support comprising the steps of: 

5 Identifying a plurality of phase states from a respiration cyde; 

measuring airflow from a subject during the sutJiecfs respiratory cycle; 
calculating phase-state probabilities for said phase slates with a measure of 
said airflow^ and 

determining a phase in the subject's respiratory cycle from the calculated 
10 phase-state probabilities of said phase states. 

2. The method of claim 1 wherein said calculating step Includes the sut>«teps 

of: 

calculating prior probabilities for said phase states; and 
15 calculating observed probabilities for said phase states given said measure of 

airflow. 

3. The method of claim 2 wherein said prior probabilities are functions of a 
phase state duration relative to respiratory cycle duration. 

20 

4. The method of claim 3 wherein said prior probabilities are a further function 
of previously calculated phase-state probabilities. 

5. The method of claim 3 wherein said observed probabilities are defined by a 
25 comparing function that relates said phase states to sub-ranges selected from a 

range of measures of an airflow cycle. 

6. The method of claim 5 wherein said comparing function changes as a 
result of changing airflow cycles of the sul^ecL 



30 



7. Ihe method of claim 3 wherein said detemiining a phase step involves 
comparing said calculated phase-slate probabilities to identify a calculated phase- 
state probability that Is more probable than other calculated phase-state probabilities. 
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8. The method of claim 3 wherein said identifying step involves partitioning a 
respiration cycle in relation to a peak flow. 

9. The method of claim 3 wherein said identifying step is repeated to adjust 
5 said phase states to account for the subject's changing respiratory characteristics as 

assessed by a measure of said airflow over time. 

1 0. The method of claim 3 wherein said step of calculating phase-state 
probabilities Involves an evaluation of the prior prdoabilities and the observed 

10 probabilities as a function of Bayes' theorem. , 

1 1. The method of daim 3 further comprising the step of adjusting ventilatory 
support based upon resulte from said determining a phase step. 

15 12. The mettiod of daim 3 wherein said step of calculating phase-state 

probabilities includes adjusting a phase-state probability based upon state 
succession. 

1 3. The method of daim 3 wherein said step of determining a phase involves 
20 tiie sub-step of calculating a continuous phase variable from said phase-state 

probabilities. 

1 4. The method of daim 1 wherein said step of measuring airflow indudes 
the sub-step of smoothing said airflow by performing a low pass filtering operation on 

25 said airflow. 

15. An apparatus for determining a subject's respiratory phase during 
ventilatory support comprising: 

a ventilator to provide a controlled supply of breattiable gas to a subject's 

30 ainway; 

a pressure transducer for generating an air flow signal representative of 
respiratory air flow from the subject, and 

a processor witii programmed instructions for conti-olling tiie steps of: 
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identifying a plurality of phase states of a respiration cycle; 

determining a measure of air flow from said air flow signal; 

calculating phase-state probabilities for said phase states with the 
measure of air flow; and 
5 detennining a phase in the subject's respiratory cycle from the 

calculated phase-state probabilities of said phase states. 

16. The apparatus of claim 15 wherein said calculating step Includes the sub- 
steps oh 

10 calculating prior probabilities for said phase slates; and 

calculating obsen/ed probabilities for said phase states given said measure of 

airflow. 

17. The apparatus of claim 16 wherein said prior probabilities are functions of 
15 a phase state duration relative to respiratory cycle duration. 

18. The apparatus of claim 17 wherein said prior probabilities are a further 
function of previously calculated phase-state probabilities. 

20 19. The apparatus of claim 17 wherein said observed probabilities are defined 

by a comparing function that relates said phase states to sub-ranges selected from a 
range of measures of an airflow cycle. 

20. The apparatus of claim 19 wherein said comparing function changes as a 
25 result of changing airflow cycles of the subject. 

21 . The apparatus of claim 17 wherein said detennining a phase step 
involves comparing said calculated phase-state probabilities to identify a calculated 
phase-state probability that is more probable than other calculated phase-state 

30 probabilities. 

22. The apparatus of claim 17 wherein said identifying step involves 
partitioning a resplrafion cycle in relation to a peak flow. 
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23. The apparatus of claim 17 wherein said identifying step is repeated to 
adjust said phase states to account for the subject's changing respiratory 
characteristics as assessed by a measure of said airflow over time. 

5 24. The apparatus of claim 1 7 wherein said step of calculating phase-slate 

probabilities involves an evaluation of the prior probabilities and the obsen^ed 
probabilities as a function of Bayes' theorem. 

25. The apparatus of claim 17 further comprising instaictions for controlling 
1 0 an adjustment to ventilatory support based upon results from said detemiining a 

phase step. 

26. The method of claim 17 wherein said step of calculating phase-state 
probabilities includes adjusting a phase-state prot^ability l>ased upon state 

IS succession. 

27. The method of clam 1 7 wherein sa/d step of detenmining a phase 
involves tiie suthstep of calculating a continuous phase variable from said phase- 
state probabilities. 

20 

28. The metiiod of claim 15 wherein said step of determining a measure of 
airflow includes the sub-step of smoothing data from said airflow by performing a low 
pass filtering operation on said data from said airflow. 

25 29. An apparatus for determining a subject's respiratory phase during 

ventilatory support comprising: 

a means for providing a controlled supply of breathable gas to a subjecTs 

ainway; 

a means for generating an airflow signal representative of respiratory air flow 
30 from tiie patient; and 

a control means with programmed instructions for controlling tiie steps of. 
identifying a plurality of phase states of a respiration cycle; 
detennining a measure of air flow from said air flow signal; 
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calculating phase-state probabilities for said phase states with the 
measure of air flow; and 

determining a phase in the subject's respiratory cycle from the 
calculated phase-state probabilities of said phase states. 

5 

30. The apparatus of claim 29 wherein said calculating step includes the sub- 
steps of : 

calculating prior probabilities for said phase states; and 

calculating observed probabilities for said phase states given said measure of 

10 airflow. 

31 . The apparatus of claim 30 wherein said prior probabilities are functions of 
a phase stete duration relative to respiratory cycle duration. 

15 32. The apparatus of claim 31 wherein said prior probabilities are a furtiier 

function of previously calculated phase-state probabilities. 

33. The apparatus of claim 31 wherein said observed probabilities are defined 
by a comparing function that relates said phase states to sub-ranges selected from a 

20 range of measures of an airflow cycle. 

34. The apparatus of claim 33 wherein said comparing function changes as a 
result of changing airflow cycles of the subject 

25 35. The apparatus of claim 31 wherein said detennining a phase step 

involves comparing said calculated phase-state probabilities to identify a calculated 
phase-state probability that is more probable tiian other calculated phase-state 
probabilities. 

30 36. The apparatus of claim 31 wherein said identifying step involves 

partitioning a respiration cycle in relation to a peak flow. 
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37. The apparatus of claim 31 wherein said identifying step is repeated to 
adjust said phase states to account for the subject's changing respiratory 
characteristics as assessed by a measure of said airflow over time. 

5 38. The apparatus of claim 31 wherein said step of calculating phase-state 

probabilities involves an evaluation of the prior probabilities and the observed 
probabilities as a function of Bayes' theorem. 

39. The apparatus of claim 31 further comprising instructions for controlling 
10 an adjustment to ventilatory support based upon resulte from said detenmining a 

phase step. 

40. The method of claim 31 wherein said step of calculating phase-state 
probabilities includes adjusting a phase-state probability based upon state 

13 succession. 

41 . The method of claim 31 wherein said step of detenmining a phase 
involves the sub-step of calculating a continuous phase variable from said phase- 
state probabilities. 

20 

42. The method of claim 29 wherein said step of determining a measure of 
airflow includes the sub-step of smoothing data from said airflow by perfonning a low 
pass filtering operation on said data from said airflow 
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